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The generalized parton distributions are non-perturbative objects, which encode information 
on long distance dynamics in a number of exclusive processes. They are hybrids of conventional 
parton densities, distribution amplitudes and hadron form factors. We give a brief review 
of theoretical developments in understanding of their properties, higher order perturbative 
effects, power corrections, and experimental observables where they are accessible. 



1 Exclusive QCD: from form factors to generalized parton distributions 

The main issue of hadronic physics is the determination of the hadron's wave function, which 
encodes a fundamental information on its structure and dynamics of elementary constituents 
forming the hadron. The cleanest indirect access to the former is achieved in lepton-hadron 
experiments, which do not suffer from difficulties intrinsic to hadron-hadron reactions due to 
initial and final state strong interaction. Chronologically, these are electromagnetic form factors 
extracted from the elastic lepton scattering off hadrons IN — > £'N', which provided charge and 
magnetization distributions within the nucleon. In a microscopic picture the form factor arises 
in the parametrization of an off-forward matrix element of the (local) quark electromagnetic 
current and is related to the wave functions by the Drell-Yan-West overlap formula 



F(A 2 ) = (P 2 |V>(0)A/>(0)|Pi) ~ / dx / dk ± W* (x, k ± + (1 - x)A ± ) & (x, k ± ) . (1) 



Here x is the longitudinal momentum fraction of the struck parton in the nucleon and k± / k± + 
(1 — x)A^ being its transverse momentum before/after the interaction with the probe. P stands 
for an appropriate Dirac matrix matching the quantum numbers of the external states. 

The measurements of P(A 2 ) gave a first insight into the composite substructure of hadrons, 
which was unraveled in its full depth in the pioneering deeply inelastic experiments £N — > i'X at 




SLAC. The Bjorken scaling observed in this reaction, has found an explanation in the Feynman's 
naive parton model that expresses the cross section in terms of a probability density f(x) to 
find a partons with momentum fraction x in a hadron. The emerging probabilistic picture relies 
on the fact that the constituents in high energy processes behave as a bunch of noninteracting 
quanta at small space-time separations. The rigorous field theoretical basis is built on the 
asymptotically free QCD and the use of the factorization theorems, which give the possibility to 
separate the contributions responsible for physics at large and small distances involved in any 
hard reaction. The field-theoretical content of the large distance contribution f(x) is given by 
the Fourier transform of the forward matrix element of a non-local light-cone operator (n 2 = 0) 

f{x) = J d\e- iXxP n (P\$(\n)ri;{0)\P) oc J dk x <?* {x, k±) & (x, fcj.) , (2) 

where, e.g. r ~ 7+, 7+75, Here we gave also an overlap type representation in terms of 

hadron wave functions, where the lowest Fock components have to be kept. 

A more direct access to the wave function of a hadron (e.g. a meson) is provided by exclusive 
production of the latter in the final state, e.g. 7*7* — > vr°. At large momenta of the virtual 7- 
quanta one probes an integral of the lowest Fock component of \P, namely 

iXxP-n 



dk ± ${x,k ± ) = J d\e- lAxy - n (P\ijj{Xn)r^(0)\0) + •••. (3) 

Recently, we have witnessed the progress in the identification of new non-perturbative charac- 
teristics that contain exhaustive information on the nucleon wave function. This was triggered 
by finding a new class of parton distributions, the so-called generalized parton distributions 
(GPDs), which interpolate between the non-perturbative functions we have discussed so far. 
Similarly, to all previous characteristics, they are defined as Fourier transforms of non-local 
quark/gluon operators sandwiched between the states with different momenta, e.g. for quark 
fields, 

A(x,r 1 ,A 2 ) = f d\e- iXx{Pl+P2) - n (P 2 \^(Xn)ri;(-Xn)\P 1 }, (4) 



and similarly for gluons. Altogether, at leading twist level there are three quark and three gluon 
twist-two operators. The parametrization of e.g. the non-local vector current ip^j^ip involves two 
leading twist functions H and E analogous to Dirac and Pauli form factors. Similarly, in the 
parity odd sector we have H and E, which correspond to the axial-vector and pseudoscalar form 
factors, etc. While being known for a while EJ, GPDs have attracted an essential attention only 
recently, after it was realized that they can shed some light onto the spin content of the nucleon. 
Namely, the second moment ( Ji's sum rule) of H + E is related 13 to the total orbital angular 
momentum fraction, J, carried by partons in the nucleon, lim^^o / dxx(H + E) = ^J. 

The cleanest reaction that gives access to GPDs is the deeply virtual Compton scattering 
(DVCS), 7* (q + f ) N -> 7 (q - f ) N ( E ^ A ) • In addition to this, they enter as a soft 

function in a variety of hard meson production, 7* ./V — > MN' , and diffractive processes. In QCD 
the Fourier transform (Ft) of the hadronic part of these processes at large momentum transfer 
(the answer to 'How large is it really?' heavily depends on the underlying strong interaction 
dynamics and presently this can mostly be judged from experiment only) is decomposed as 

(N'\T{j^z)j u (0)}\N) ^T^ u j\x C(x,OA(x,r ] )+T^j\x C 3 (x, £)A 3 (x, 77) + 0(<T 2 ),(5) 

{N'M\j^(z)\N) ^T^j\xj\y <f)(y)a s C(y, x, £)A(x, rj) + 0(Q~ l ), (6) 

for DVCS0'§! and exclusive meson production! amplitudes, respectively. Here C stands for the 
hard scattering subprocess, while A and cj) stand for GPDs and meson distribution amplitudes 




(DAs), respectively. Here T is a Lorentz tensor. The scaling variables are £ = —q 2 /q ■ P as well 
as 77 = g • A/g • P (~ — £), and Q 2 = —(q + A/2) 2 . In Eq. (||) we have kept the power suppressed 
contributions, T 3 V oc A±/Q, from twist-three GPDs A 3 . 

2 Properties and models for GPDs 

In different region of the phase space GPDs share common properties with conventional parton 
densities for |x| > 77, and distribution amplitudes for \x\ < 77. GPDs can not be interpreted as 
densities in general but rather as interference terms between wave functions of incoming and 
outgoing hadronsB. From the operator definition (|J) it follows that the j th moment of GPDs 
A = {P, E, E}, Aj(r]) = J dxx^~ 1 A(x, 77), is a polynomial of order j in skewedness 77 while for 
P is of (j — l) st only. The sum H -\-E also obeys the latter property so that the rj independence 
of the Ji's sum ruleB is a particular example. To preserve the polynomiality condition of GPDs 
the parametrization in terms of spectral, or double distribution (DD), function reads 

x, 77, A 2 ) = J ^dyj dz6(y + r,z - x) j X ~ j (y, z, A 2 ). (7) 

While the last definition coincides with the original one introduced by0i, the former ill differs 
from it and respects the polynomiality condition alluded to above. An alternative solution given 
in Ref. □ consists of keeping for A the representation in terms of DDs of the second line and 
adding an independent function D{x/rf) to it so that this term produces the missing r/^-term in 
the j th moment. The inverse transformation of DDs in terms of GPDs has been derived in Ref. 
Ill and was found in ^1 to be in one-to-one correspondence with the Radon transformation. 

Let us consider a particular example of GPDs. One assumes a factorizable ansatz of (x, rj) 
and A 2 dependence: H(x,rj, A 2 ) = ElA 2 )H(x,r]), with H(x,rt) expressed in terms of a double 
distribution F(y,z) that is modeled" as 

F(y, z) = {Aq(y)9(y) - Aq{-y)6{-y)} *(\y\,z), n(y, z) = (8) 

and obeys the F(y,—z) = F(y,z) symmetry 0. For illustration purposes let us discuss the 
isotriplet combination. Assuming an SU(2) symmetric sea Aq^(y) = Au(y) — Ad{y) = 0, 
we have Aq^(y) = Au va \(y) — Ad va x(y) and model it by a semi-realistic ansatz Aq( 3 \y) = 

9 A r(i)T(i-n) ^ an< ^ 9 A = 1-26. This allows to get a simple analytical representation 
H^(x, V ,A 2 ) _ (1-f) ^{x + V 



9(x > —77) ^ y~~^ J {j] 2 — x + (2 — 71)77(1 — x)^j — (77 — ► —77) 



5A P(A 2 ) 7/3 I v ' n \l + 7] 

'(9) 

where a dipole parametrization P(A 2 ) = (1 — A 2 /m 2 ^) 2 with m 2 A = 0.9 GeV 2 is used. One sees 
that H^'(x. rj) vanishes for x < —77. Other models arise from computations in the framework 
of the bagl!3 and the chiral quark soliton model §0 or based on the overlap representation B. 



3 Power suppressed corrections 

The leading twist-two predictions (|5|) (with A 3 = 0) and (J6|) are affected by a number of cor- 
rections, with power suppressed effects in Q being one of them. They have been addressed in 
detail recently for the generalized Compton amplitude. In Ref.lll, using a QCD improved parton 
picture of Ellis et al., it was shown for the example of a scalar target that the twist-three con- 
tributions induce gauge restoring pieces (lacking in the twist-two approximation) in theLorentz 
structure of the Compton amplitude. A complete OPE-based analysis has been done in&3 where 



generalized Wandzura-Wilczek relations have been found. The twist-three GPDs are expressed 
by these relations in terms of twist-two ones as well as i nteraction dependent antiquark-gluon- 
quark correlations, A 3 = W (g> A + (ipGip). In Refs. I 17 H 18 I it was observed that for the DVCS 
kinematics, r\ ~ — £, the Wandzura-Wilczek relations alluded to above develop a singularity, 
however, the latter does not show up in the cross section in which the twist-three functions en- 
ter in a specific combination. In an earlier attempt El a generalized Wandzura-Wilczek relation 
has been derived by neglecting all twist-three operators, and is, therefore, erroneous. 

Since the cross sections for the reactions in question are peaked at low momentum transfer, 
mass corrections may play a vital role in the confrontation of theoretical predictions with ex- 
perimental measurements. As a first step towards this direction we have resummedB the target 
mass corrections, Ylk c k {M 2 /Q 2 ) k , stemming from the trace terms of the twist-two operators. 
The same can be done in the twist-three sector and will provide effects of the same order. The 
twist-four sector awaits its unraveling since the renormalon analysis demonstrated a potentially 
sizable effect coming from the latter already for the moderate momentum transferee. 

4 Perturbative corrections 

Leading order calculations in QCD are known to be unreliable and they are strongly affected by 
perturbative corrections. For the cases at hand, the coefficient function is a series in coupling 

C = C + + • • • . (10) 

The one-loop correction C\ to th e han dbag approximation of the generalized Compton scattering 
amplitudes have been derived inrT 2 o. For the hard meson production with quark dominated 
parton subprocess, C\ has been recently extracted in Ref.El making use of the known result for 
the pion form factor. The analyses have shown that they can produce a very large modification 
of the leading order (LO) predictions @@. 

Yet another source of corrections stems from the evolution of GPDs (and DAs), which is 
governed (at twist-two level) by the equation (with r\ = 1 for DAs) 

d f 1 a / a \ 2 

A(x,r,)= J dyK(x,y,ri)A{y,Ti), with K = ^-K + K x + . . . . (11) 



din Q' 

The one loop kernels for all channels have been computed by Lipatov et al. in Ref. ^ and 
recalculated by a number of groups in the last few years of XX th century (see i for a complete 
list of references). The two- loop Quantities have become available first in the form of local 
anomalous dimensions from Refs. E3't3'E3 and later have beenconverted into the momentum 
fraction form in 123 by means of an extensive use of conformalIHa and supersymmetric Eil Ward 
identities. Note that the flavour nonsinglet kernel has been known for quite a while fromilil'ii. 

The evolution effects have been studied by methods based on the orthogonal polynomial 
reconstruction 00, extended inS to NLO accuracy, and the direct numerical integration EL 
The resulting two-loop effects have been found to alter the LO evolution in a percentage range. 



5 Observables 

In the electroproduction process of a real photon the DVCS signal is strongly contaminated by 
the Bethe-Heitler (BH) process, a oc |DVCS + BH| 2 . However, individual terms in this sum 
are distinguished by peculiar dependencies on the lepton charge, lepton A and hadron spin, and 
the azimuthal angle ip of the outgoing photon. The best quantity to access the GPDs is, of 
course, the DVCS-BH interference since the GPDs enter linearly in it. Its precise experimental 



extraction is only possible on machines having lepton beams of both charges, a + — a~ oc DVCS- 
BH. This allows a clean separation of Fourier components w.r.t. azimuthal angle, DVCS-BH 
oc ^ m [c m cos(m^) + \s m sm(mip)]. For a lepton beam with a given charge, other observables 
(like diverse spin asymmetries, which isolate the DVCS-BH interference at leading power in 
1/ Q-expansion) will be contaminated by (sizable) power suppressed (higher twist) effects. For 
instance, in the cross section a, the coefficient of cos / sin(<^) stemming from the twist-two DVCS- 
BH will get corrected by |DVCS| 2 at twist-three level, etc. Note however, that the twist-three 
effects do not induce the same angular dependence in DVCS-BH and |DVCS| 2 , separately, which 
is already generated by leadingtwist contributions, i.e. in DVCS-BH they produce cos / sin(2y?) 
only (for a numerical study seeE3), so that the coefficient of cos / sin(^j) will be modified at twist- 
four level 0. The experimental isolation of different Fourier components can be done by forming 
appropriate azimuthal asymmetries or weighting data with corresponding angular functions 
cos / sm(m<p). As distinguished from DIS, in DVCS one can access tensor iduons, which are 
not contaminated by quarks at twist-two. They show up at one-loop order Ej-Ej and generate 
a specific cos / sin(3(p) azimuthal angle dependence in the cross section Provided we have 
separated DVCS-BH interference alone, the single spin asymmetries make it possible to extract 
the imaginary part of the DVCS amplitude and thus give acc ess to the shape (at LO in a s in 
complete analogy to DIS) of GPDs on the diagonal x = £, e.g.cirE 2 ! 

da^-da^ m y 2 (2 - s B )(2 - y) f x B „ ^ A 2 1 

|Aj_| smy>Im< F{H + &m n — 9 F 2 £ > , 



dx B dQ 2 d\A 2 \dip tt ^T^Q 5 A 2 1 X| r [ 2 - x B ±M 2 ' 

'(12) 

where A = C®A, see Eq. (|5|), Fi, F2 and Gm = F1+F2 are Dirac, Pauli and magnetic Sachs form 
factors, and the transverse momentum squared is A^ = ((1 — xb)A 2 + x B 2 M 2 ) /(l — x B /2) 2 . 

In hard exclusive meson production one measures products of off- forward amplitudes, e.g. 
the asymmetry of the 7r + -production off a transversely polarized proton reads EJea 

da ± -da l _ 4tt/ 2 4 \A ± \ ^^ftfa)* c(3)) (13) 
d\A 2 \d^ ~ m 9 Q 6 2 — xb M sm ^ Lm \ n t /' W 

with (Tt,£) = <f^£&JJLE), see Eq. Cross sections for production of other meson species 
can be found inEJ'c3ac3c§cZI. Note that quark transversity GPD does not show jup in the meson 
production due to preservation of the chiral symmetry in perturbation theory The current 
experimental situation is discussed inc3 andH for HERMES and JLab settings, respectively. 
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